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Using the structural data of the La^CnO^ compound both in the LTT phase and in the isotropic 
phase we have derived an effective t — J model with hoppings t and superexchange interactions 
J extended up to fourth and second neareast neighbors respectively. By numerically studying this 
hamiltonian we have then reproduced the main experimental features of this HTc compound: d-wave 
superconductivity is stabilized at small but finite doping 5 > 6% away from the antiferromagnetic 
region and some evidence of dynamical stripes is found at commensurate filling 1/8. 

PACS numbers: 74.20.Mn, 71.10.Fd, 71.10.Pm, 71.27. +a 



The microscopic mechanism of high-temperature su- 
perconductivity (HTc) is probably the most important, 
but still open problem in condensed matter physics. Af- 
ter the recent discovery of MgBi [EJ - a phonon like 
superconductor with Tc ~ 39K, and the fullerene su- 
perconductivity temperature 0] has been substantially 
enhanced in a compound where certainly the electron- 
phonon coupling cannot be neglected, there have been 
increasing expectations ]3j that the conventional type of 
electron-phonon BCS superconductivity may explain all 
HTc materials || . On the other hand, from the theoret- 
ical point of view, a wide range of numerical techniques 
have recently led [Q] to rather consistent evidence that 
strongly correlated models, such as Hubbard like ||, and 
t - J models @ § may show d-wave superconductivity 
at low temperature, evidence that is rather remarkable 
because no electron-phonon or explicit attractions are in- 
cluded in these models. 

In order to contribute to the solution of these contro- 
versial aspects between theories and experiments and to 
understand the role of the strong electron correlation, in 
the present work we try to reduce the distance between 
the physics of rather abstract model hamiltonians and a 
consistent ab-initio derivation of experimental properties 
of HTc material. We consider the most popular com- 
pound for HTc superconductivity - LaiCuO^. We de- 
rive, and then study by Quantum Monte Carlo methods, 
an effective model in which the electron-phonon interac- 
tion is neglected. Our purpose is to understand whether 
the physics of an effective model hamiltonian may ex- 
plain properties of real materials that are not understood 
with conventional methods of band theory and/or BCS 
mechanism of superconductivity. To our knowledge our 
work represents the first attempt to explain the physics 
of HTc superconductors, starting from first principle cal- 
culations. 

The LaiCuO^ compound is known to be supercon- 
ductor in a range of doping between 6% and 30%, and 
by interstitial substitution of Nd (in place of La) at 



commensurate filling 1/8, static incommensurate peaks 
in the magnetic structure factor were observed by neu- 
tron scattering experiments H and, correspondingly, 
the superconducting transition temperature was strongly 
suppressed. The Nd-substitution enhances the spatial 
anisotropy and is believed to favor the formation of 
stripes: one dimensional hole-rich patterns in the CuO 
planes, separating antiferromagnetic (AF) domains with 
opposite direction of the AF order parameter, namely 
the next nearest neighbor Cu-spins across the stripe are 
antiparallcl. This feature easily explains the mentioned 
neutron scattering incommensurate peaks, within the as- 
sumption that the stripes arc half filled (half an hole per 
Cu site density along the stripe). In fact this implies that, 
at a given hole doping S, the incommensurate peak in the 
magnetic structure factor S(q) shifts by 2nd from the 
AF wavevector (tt, tt) — > (tt, tt — 2tt5) (for equally spaced 
stripes parallel to the x-axis direction), correspondingly 
the charge structure factor N(q) shows up an incommen- 
surate peak close to the T point for q = (0, 4ir5) ||. 

We consider the extended t — J one-band model on a 
finite lattice with N sites: 



R,fl,<7 



where c. 



(1 — rii t ff), and is the electron spin 



operator on site i, whereas the sum run over the lat- 
tice sites R and corresponding neighbors R + r M , deter- 
mined by vectors r M shown in Table I. Previous work 
have shown the key relevance of hoppings extended be- 
yond nearest neighbors to reproduce the phenomenology 
of the Cuprates . In this table we microscopically de- 
rive the values of the effective parameters t M and J M de- 
termined from the results of electronic structure calcu- 
lations by standard LDA (TBLMTO [jlOj) and LDA+U 
method for the anisotropic case (with Nd substi- 
tution) and the high-symmetry tetragonal case, using 
known structural data |T^ | . The effective hopping param- 
eters tfj, were calculated by the standard least-square fit 
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TABLE I: Values of the effective parameters in eV. 





Anisotropic 




Isotropic 




T 






J* 




(1,0) 


0.105 


0.425 


0.109 


0.486 


(0,1) 


0.111 


0.466 


0.109 


0.486 


(1,1) 


0.016 


0.014 


0.016 


-0.086 


(2,0) 





0.036 





-0.006 


(0,2) 





-0.064 





-0.006 


(2,1) 





-0.001 








(1,2) 





0.046 









procedure to the bands obtained in LMTO calculations. 
The effective exchange parameters J p were calculated us- 
ing the formula derived by A. Lichtenstein Jl3|, [lj] , where 
the second derivative of the total energy as a function of 
the value of the angle between spin directions on two 
atoms is calculated via eigenvalues and eigenfunctions 
obtained in electronic structure calculations. Those elec- 
tronic structure calculations were done using LDA+U 
method which allows to obtain antiferromagnetic insulat- 
ing ground state for the undoped cuprate. By contrast, 
the standard LDA method leads to a nonmagnetic metal- 
which is obviously inconsistent with 
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experiments. The Coulomb parameters U and J used in 
LDA+U calculations were obtained in constrained LDA 
calculation || (U=ll eV, J=l eV). 

The Nd-substitution for La in La,2CuOi results in 
anisotropic Low Temperature Tetragonal (LTT) struc- 
ture JT2| (Fig.|l|), where the CuOq tetrahedra are tilted 
via rotation around [1,0] axis (and [0,1] axis in the next 
plane, so that the total crystal structure preserves tetrag- 
onal symmetry). As a result of this tilting there are two 
kinds of oxygen atoms, 01 which form 180 degree Cu- 
O-Cu bonds in [1,0] direction and are in the Cu atoms 
plane, and 02 atoms, which are shifted out of the plane 
and form Cu-O-Cu bonds in [0,1] direction with the de- 
creased bond angle. Remarkably the strongest nearest 
neighbor hopping is in the direction [0,1] where the oxy- 
gen atoms are out of the Cu-Cu plane (see Fig. [l]). This 
seems counterintuitive, because for this direction the Cu- 
O-Cu bond angle is less than 180 degree and hence the 
effective d-d hopping must be decreased comparing with 
the straight Cu-O-Cu bond in [1,0] direction. Solution 
of this puzzle is the difference in electrostatic potential 
on 01 and 02 crystallographic positions which results in 
the energy of 2p-orbitals of out-of-plane 02 atoms being 
1.05 eV higher than for the in-plane 01 atom. As the 
effective d-d hopping is proportional to t^ d /A p d, where 
t p d is effective p-d hopping parameter and A p d is energy 
separation between Cu-3d and 0-2p orbitals, then A pc i is 
significantly smaller for out-of-plane 02 oxygen resulting 
in a larger value of the effective d-d hopping parameter. 
Another effect of anisotropic crystal structure distortion 
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FIG. 1: Structural data for the HTc materials considered 
(see text). The displacements of the ions in the anisotropic 
LTT phase is indicated by the arrows. 



is the long-range nature of model hamiltonian hopping 
parameters. For the isotropic case only two shells of 
nearest neighbors were needed and the value of for the 
third neighbor is negligible. Instead for the anisotropic 
case the hopping parameters have to extend at least up 
to fourth neighbors for a satisfactory fit of the LMTO 
bands. We must note also, that the usual scaling of the 
effective exchange parameters as is not valid here. 

In order to study this model hamiltonian we use the 
recent QMC methods, that have been shown to work 
very well for the standard t — J model jL6|, allowing to 
obtain very accurate values of energies and correlation 
functions for a wide range of J/t values. The initial vari- 
ational wavefunction is chosen to be the most general 
Jastrow-BCS real wavefunction i/jq = J\BCS > pro- 
jected onto the subspace with no-doubly occupied sites 
and fixed number of particles. After the latter projection 
the wavefunction is totally symmetric, and a large num- 
ber of variational parameters can be optimized consis- 
tently with the symmetries of the model. This is the first 
step p = of our procedure (SR) |l6|, [lTj] (e.g. 46 varia- 
tional parameters for the 8x8 anisotropic case) , which is 
in principle convergent for large number p of iterations. 
The scheme is based on the fast convergence properties 
of the Lanczos technique, that allows to improve remark- 
ably the accuracy of the best variational guess with few 
iterations. Starting from the p— Lanczos step wavefunc- 
tion it is also possible to further improve the variational 
energy and especially the quality of the correlation func- 
tions, by a further correction scheme, such as the fixed 
node (FN) |l8| or the more accurate stochastic reconfig- 
uration method (SR) [Tij ]. A comparison of the quality 
of our approximation, as compared with the conventional 
Fixed node technique |L8| is shown in Tab. [n]. 

We concentrate our study in the low-doping region 
and especially at the important filling S — 1/8 where 
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TABLE II: Variational energy (eV) per Cu site for the 
anisotropic t — J (Tab. |j|) model for the best variational FN 
and SR (p = 2) techniques applied to the p = 1 Lanczos wave- 
function. The a = are estimates of the exact zero variance 
fl6l. Error bars are in brackets. 



energies 



Anisotropic 



Isotropic 



# Holes 


VMC 


FN 


SR 


(7 = 





-0.06446(1) 


-0.06614(4) 


-0.06637(4) 


-0.06665(17) 


4 


-0.12717(3) 


-0.1368(1) 


-0.1379(1) 


-0.1402(5) 


6 


-0.15623(4) 


-0.1679(1) 


-0.1701(2) 


-0.1744(14) 


8 


-0.18336(4) 


-0.1972(1) 


-0.1988(3) 


-0.2025(7) 


10 


-0.20859(4) 


-0.2237(1) 


-0.2265(3) 


-0.2316(6) 


16 


-0.27446(6) 


-0.2929(1) 


-0.2956(4) 


-0.3053(11) 


26 


-0.35030(6) 


-0.3684(1) 


-0.3710(4) 


-0.3780(12) 



static stripes were observed. As shown in Fig. [2] there is 
a clear evidence, especially in the anisotropic case, of 
an incommensurate peak that was not present at the 
variational level, and is sharpening as the accuracy of 
the calculation approaches the low energy limit. Cor- 
respondingly the charge structure factor N(q) does not 
seem to be enhanced in this limit, but only a cusp sin- 
gularity at the Tranquada's wavevector appears consis- 
tent with our data. The strongest anisotropy appears in 
the nearest neighbor hoppings (see Tab. 1) and, corre- 
spondingly the peak in the magnetic structure factor, is 
shifted by (0, 2nS), consistent with the Tranquada's pre- 
dictions. Within the "stripe picture" this is compatible 
with stripes along the [1,0] axis, a direction where the 
nearest neighbor hopping is much less favored compared 
with the other direction. Thus, in order to optimize the 
kinetic energy in the y-direction, charge-fluctuations per- 
pendicular to the stripes are required. These fluctuations 
suppress any static response in the charge structure fac- 
tor, but leave a sizable effect in the spin provided the 
spins across the hole are antiparallel [|l9] , as we have ver- 
ified in this case. In the following we name this feature 
the " dynamical stripe" , a genuine effect of strong corre- 
lation. 

As shown in Fig. ^|, the anisotropy clearly enhances the 
dynamical stripe fluctuations, since, in the anisotropic 
case, a more resolved incommensurate peak for S(q) is 
found. Some small effect is seen also in the isotropic case, 
consistent with the DMRG findings jjj that t' jt < may 
stabilize stripes. Probably at very low energy (larger 
p) evidence of true static stripes can also be found at 
this particular filling. For all other dopings, though in- 
commensurate magnetic peaks are still present, they are 
much less sharp, and the value of S(q) at the maximum 
is much below the corresponding p = variational value. 

Regarding superconductivity, we have studied the 
pairing correlation functions A^' 1, ' (r) = (SiJ rr ^S i J) . Here 

Sj ^= c[ tc| + ^ — c| , c| + . creates an electron singlet pair 
in the neighboring sites (i,i+fi). Off-diagonal long-range 
order is implied if Pd = 2 lim,. 




FIG. 2: Spin and charge structure factors for the cases of 
Tab. I vs. improved accuracy in energy: VMC, SR p=l, SR 
p=2, respectively. T = (0,0). X = {tt,tv), M = (tt,0), M = 
(0, 7r) and the arrows indicate the Tranquada's wavevectors 
(see text). 
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FIG. 3: Pd for various dopings and improved accuracy in 
energy, VMC, SR p=l, and SR p=2, respectively, for the 
cases of Tab. |. Lower panels: the value of Pd was obtained by 
d— wave-averaging the pairing correlations at the maximum 
distance ~ 4\/2. Upper ones: Pd refers to the maximum 
distance (= 4) along the x,y directions. 



nite in the thermodynamic limit. 

For the 8x8 system, that we have studied for several 
dopings (see Fig. || lower panels) size effects are accept- 
able, at least close to optimal doping as shown in |^| for 
the t — J model. The existence of a finite Pd and corre- 
spondingly the absence of antiferromagnetic long range 
order is an experimental fact for S > 6%. Both features 
are remarkably well reproduced by our calculations (see 

ht Fig. |). In particular, at half filling, the magnetic 
structure factor dramatically increases with respect to 
the spin liquid variational wavefunction J\BCS), con- 
firming the existence of AF long range order, as widely 
accepted for the weakly frustrated Heisenberg model. 
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FIG. 4: Right panel: S(n, n), as a function of doping 8 for the 
8x8 anisotropic (Tab. |) t — J model. Left panels: S(tt, ty) and 
Pd for larger size for the standard t — J model (including the 
nearest neighbor density-density interaction, which however is 
irrelevant at small doping) and the anisotropic case indicated 
by the arrow. In the latter case the J/t value is defined by 
the average nearest neighbor ratio. 



ing completely the electron phonon coupling, within an 
ab-initio derivation of the effective superexchange and 
long-range hopping couplings for an extended t— J model. 
Several interesting features comes out from our calcula- 
tion that deserve experimental confirmation: the finite 
critical doping S c required to stabilize superconductivity, 
appears a band structure effect. In principle antiferro- 
magnetism and superconductivity may coexist at small 
enough doping, provided the long range hoppings are sup- 
pressed. The anisotropy and long range hoppings appear 
to be compatible with a sizable superconducting order 
parameter, as long as the stripes remain dynamical [ p0[ . 
Static stripes are not clearly stabilized in a model in 
which the electron-phonon interaction is disregarded. 

This work was partially supported by Russian foun- 
dation for Basic Research grant RFFI-01-02-17063 and 
by MIUR-COFIN2001. We acknowledge T.M. Rice 
for stimulating the present scientific collaboration and 
E. Dagotto for useful discussions. 



Conversely, as soon as the doping is tinily increased, the 
AF magnetic structure factor remains very close to the 
spin liquid variational reference, clearly indicating ab- 
sence of AF long range order. This is the main result of 
our paper and is due to the long range couplings of our 
effective model, since at small doping the t — J model 
without long range couplings is expected to be both an- 
tiferromagnetic and superconductor ||. In order to un- 
derstand this effect, we have extended the calculation to 
larger size N — 242 at small doping both for the standard 
t — J model at J/t = 0.2 and J/t = 0.4 and our extended 
t — J model, in the anisotropic case. As it is clearly seen 
in the left panels of Fig. [|, the role of the long range hop- 
pings is crucial to destroy both antiferromagnetism and 
superconductivity. 

Another effect is that, within our model, the anisotropy 
does not suppress the superconducting order parameter 
Pd, (see lower Fig. ^). This is rather in contrast with 
experiments, where a drastic change of T c was found in 
the anisotropic case. However in this case the pairing 
correlations are also anisotropic (see upper Fig J^) and 
Pd vs doping, computed at the maximum [1,0] or [0,1] 
distance appears to be much suppressed in the direction 
of the stripe (notice that this direction is alternatively 
x or y in neighboring CuO planes). It is possible that 
by including the clcctron-phonon coupling one can ob- 
tain further agreement with experiments. In fact we ex- 
pect that, whenever dynamical stripes are clearly formed, 
the electron-phonon coupling can considerably enhance 
the anisotropy and lead to the formation of true static 
stripes, with negligible Pd- 

In conclusion we can reproduce many aspects of the 
low-doping phase diagram of HTc compounds by neglect- 
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